J-T Chen 5-8-3 



CODE ASSIGNMENT IN A CDMA WIRELESS SYSTEM 



This application claims the benefit of U.S. Provisional Application Serial No. 
60/137462, filed June 4, 1999. 

5 

Field of the Invention 

The present invention relates generally to code division multiple access (CDMA) 
wireless communication systems. More particularly, the present invention relates to 
code assignments in CDMA wireless communication systems. 

10 

Background of the Invention 

Wireless communication systems in which wireless terminals (e.g. fixed or 
1^ mobile) conmiunicate with fixed base stations via a wireless channel are well known in 

W the art. Since the wireless spectrum bandwidth is limited, various techniques have been 

ijj 1 5 developed to allow multiple wireless terminals to share the limited bandwidth. One 
IS technique is time division multiple access (TDMA) in which wireless terminals share a 

□ single frequency but are separated in the time domain with each v^reless terminal being 

Q assigned a time slot. Another technique is frequency division multiple access (FDMA) in 

which wireless terminals are separated in the frequency domain, with each wireless 
y 20 terminal being assigned a particular frequency. 

S The present invention is concerned with another technique for sharing limited 

spectrum, namely code division multiple access (CDMA). In a CDMA system, the 
signals of all wireless terminals share a frequency band and are sent together in time. 
The signals of different wireless terminals are differentiated using codes. Each wireless 
25 terminal is assigned a unique code (S) (also known as a spreading code or signature). 
The signal (d) transmitted from a particular wireless terminal to the base station is 
combined with the wireless terminal's unique code iS to generate an encoded signal over a 
wide band of frequencies. If there are K wireless terminals (1, 2, ... AT), the signal of the 

A:^ wireless terminal is represented as {<ijt(*)} the code of the A:* wireless terminal is 
30 represented as . The transmitted signal for the li^^ wireless terminal at the /th time 

instant is represented as df^{i)Sf^ . Thus, the combined signal for all the wireless terminals 
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received at the base station is (i)S^ +• • '+df.(i)Sj^ . The codes are ideally designed to be 
orthogonal such that the signal of the kth wireless terminal may be extracted at the base 
station by multiplying the received signal with the wireless terminal's unique code Sj^ . 
Signals originating at the base station are also combined with a code prior to being 
5 transmitted to each wireless terminal. Upon receipt of the signals, the wireless terminal 
extracts the signals from the base station by multiplying its received signal by appropriate 
codes. CDMA technology is well known in the art of wireless communications and will 
not be described in further detail herein. 

One problem with CDMA systems is inter-user interference (also known as 

10 multiple access interference (MAI)) which results from the leakage of signals of one 
wireless terminal into the signal of another wireless terminal. This leakage is due to 
several factors. First, there is multipath propagation which results in the same signal 
being propagated along multiple paths with different time delays. Moreover, the codes 
are usually not perfectly orghogonal in order to be able to accommodate a large number 

1 5 of terminals. This non-perfect orthogonality itself also contributes to inter-user 

interference. Also, received signals may not all be perfectly synchronized in time, which 
also contributes to the inter-user interference. 

In order to deal with the interference problem, current systems generally perform 
complex operations at the receiver in an attempt to remove this unwanted interference. 

20 One approach of dealing with interference is single user processing, in which a receiver 
uses only knowledge of a single signal (e.g. the signal associated with a single wireless 
terminal) to process the incoming signal. A more sophisticated approach is multi-user 
processing, wherein the receiver uses knowledge of all signals in order to demodulate 
them all in a joint fashion. This multi-user appraoch is generally only performed at base 

25 stations because of the significant processing overhead which it requires. 
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Summary of the Invention 

We have realized that the interference experienced by a wireless terminal 
in a CDMA wireless communication system is a result of a synergy between the code 
assigned to the wireless terminal and the particular wireless channel through which the 
5 signals of the wireless terminal pass. Thus, a wireless terminal signal propagating 

through a particular v^reless channel will experience different interference with different 
codes. 

In accordance with the invention, codes in a CDMA wireless communication 
system are assigned based on characteristics of the wireless channels. By coordinating the 

10 assignment of codes among the v^reless terminals based on the wireless channels, the 
total interference among the wireless terminals is reduced and performance of the overall 
wireless communication system is improved. 

In accordance v^th one embodiment of the invention, an iterative code 
optimization algorithm maintains a processing set of wdreless terminals. In accordance 

1 5 v^th the algorithm, when the codes for the wireless terminals within the processing set 
are optimized, additional wdreless terminals are added to the processing set and the codes 
for the new processing set are optimized. This continues until the processing set consists 
of all the wireless terminals in the wireless system. This technique, called sequential 
packing, is advantageous because the code assignments are more likely to converge to the 

20 optimum solution when the number of vsdreless terminals is small and the optimal code 
assignment of a ^-wireless terminal system is likely to be similar to that of a (A^-1) 
wireless terminal system. 

In accordance with another embodiment of the invention, during the processing of 
each processing set a technique called target v^reless terminal isolation is used. Instead 

25 of updating all the codes for all the wireless terminals in the processing set, a target 

wireless terminal is chosen and its code is optimized prior to moving on to optimize the 
codes of the other wireless terminals in the processing set. This technique also speeds up 
convergence because before the convergence of the code optimization, most of the codes 
may have already achieved good orthogonal distribution. Disturbing all the codes would 

30 likely slow down the convergence speed. 
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In accordance with one embodiment of the optimization algorithm, the code of a 
target wireless terminal is assigned as follows. First, a random search of codes is 
performed in order to find a code which is an improvement over the currently assigned 
code. In one embodiment, the random search is performed by searching all available 
5 codes until an improved code is found. In an alternate embodiment, the random search is 
performed by searching a subset of all available codes for the best code in the subset. A 
performance analysis of these random search techniques is described in the detailed 
description. The new code foimd by the random search is then improved by performing a 
gradient search of codes. In accordance with this technique, the signal space which 
10 surrounds the new code is searched in order to further optimize the code for this wireless 
terminal. This gradient search of codes may be performed by searching the codes which 
result fi-om making a small change to the code found by the random search. The new 
code found by the random search may also be improved by performing a gradient search 
;=y of transmission delays. In accordance with this technique, the delay of transmission of 

ijj 1 5 the new code is adjusted by a small amount in order to determine the optimum delay of 
"'Z the code. 

Q The invention provides advantages over the prior art systems. Interference is 

reduced at the outset by the assignment of codes based on channel characteristics, rather 
^ than randomly assigning codes at the outset and trying to remove the interference at the 

UJ 20 receiver. This substantially simplifies receiver processing. 

^ These and other advantages of the invention will be apparent to those of ordinary 

skill in the art by reference to the following detailed description and the accompanying 
drawings. 

25 Brief Description of the Drawings 

Fig. 1 shows a system in accordance with one embodiment of the invention; 

Fig. 2 is a flowchart of the steps performed to assign codes to wireless terminals 
in accordance with one embodiment of the invention; 

Fig. 3 is a graph showing the minimum SINK distribution in an example system; 
30 Fig. 4 is a graph showing the top 2 percent performance of Fig. 3; 
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Fig. 5 is a graph showing the improvement of the probability density function of 
the minimum SINR using a first search strategy; 

Fig. 6 is a graph showing the number of search cycles required to reach certain 
performance levels; and 
5 Fig. 7 is a graph showing the improvement of the probability density function of 

the minimum SINR using a second search strategy. 

Detailed Description 

A system in accordance with one embodiment of the present invention is shown 

10 in Fig. 1. Wireless terminals llOi, IIO2, ... llOj^ are communicating with base station 
102 in accordance with the CDMA communication protocol. Broken line 104 represents 
a logical separation between forward link (i.e., base station to wireless terminals) 
components above the broken line 104 and reverse link (i.e., wireless terminals to base 
station) components below the broken line 104. Base station 102 comprises antennas 

15 IO61, IO62, ... 106/> for receiving signals from the wireless terminals 1 10. A frequency 
down converter 1 12 receives the signals from the wireless terminals 1 10 and demodulates 
the signals to extract baseband signals. Thus, the output of frequency downconverter 112 
is the combined baseband signal from all wireless terminals fc^^d^ {i)S^ +• • -^K^f^d^ (z)5^ , 

where k is the combined effect of the wireless channel of the Ath terminal and the /7th 

20 antenna element. Note that the signal model used herein does not assume intersymbol 
interference (ISI), but one skilled in the art could readily extend the model to include ISI. 
The combined signal is provided to a joint channel estimator 1 14 and a space-time 
matched filter 116. The joint chaimel estimator 1 14 receives the baseband signal and 
determines the characteristics of the v^reless channels k over which the signals have 

25 been received. 

The wireless channel is described by a multi-path propagation model. A 
mathematical description of the model which provides a basis for a further description of 
the elements of Fig. 1 v^ll now be given. In large cells with high base station antenna 
platforms, the propagation channel comprises a few dominant specular paths. In such a 

30 case, the baseband signal received by the base station 102 can be expressed as follows: 
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k=\ l=\ 

where x(/)is the received baseband signal as a function of time /; is the additive 
noise; K is the number of wireless terminals visible to the base station 102; L is the 
maximum number of propagation paths present in the system for each wireless terminal; 
5 z{o^ j ) is the steering vector of wireless terminal ^ to a signal arriving from direction ^ ; 

fif^ I is the complex time-varying amplitude of the path of the wireless terminal, 
with the path amplitude being a complex Gaussian random process including both the 
propagation loss and the signal fading caused by the Doppler spread; ?^ is the transmitted 
complex baseband signal from the h!^ wireless terminal; r^j is the propagation delay of 

10 the /'^ path of the 1^^ wireless terminal. 

In a CDMA system with the length of the spreading code A'^ equal to the number 
of chips in a symbol period, the transmitted signal 1 can be represented as the 
ill 'convolution of the data bits, the codes and the pulse-shaping ftmction. Assuming band- 

% limited signals in the system, if each continuous signal is sampled at a rate which is 

Q 1 5 higher than twice the Nyquist rate, then the convolution, correlation, and integration 

operations can be replaced by multiplication of matrices with appropriate arrangements 
of matrix elements from the sampled data. In the present embodiment, we assume P 
^ji antennas, and we oversample the signals with a factor of p (i.e., we sample the signal p 

:jQ times in each chip). In addition, there are chips in each data symbol, and Mdata 

20 symbols are considered. We also assimie that the maximum path delay plus the length of 
the pulse-shaping function is no greater than and that Q<N, where Tc is the chip 
period (if Q is much smaller than N, then the system is quasi-synchronous). After 
sampling, we can express the received signal X as 

x = Xaaga».+n (2) 

jt=i 

25 where X: P x MNp and N: P x MNp are, respectively, the received signal and the 
additive Gaussian noise; A^: P x L is given by 

A^ =[a(^A,,)a(^jfc,2)-a(^A,L)]; Bj^ilx I is given by =diag{>ff J with 



6 



J-T Chen 5-8-3 




Pk - \.Pk,\ 'Pk,2 > • • •■'Pk.Lf 5 the /,/'' element of G^.LxQp can be written as 

g(/o + - — - -'^ki) ^^1^ ^0 being the normalized pulse-shaping function. The code 
P 

matrix S,^:Qpy.{N + Q-\)p is a Toeplitz matrix with [j^ , ,^u(Qp.^)Y as its first column, 
and {Sk,x,^x^(p.^)S^,2,^x^(p.^),...,S^j^Six.(Qp-\)\ as its first row, where s^ ^,s^^,...,s^J, is the 
5 normalized code of the k!'' wireless terminal. The data matrix ^i,^•.{N + Q-\)px MNp is 
also a Toeplitz matrix with [c/^ ,^u(Np.\)^d^^_^ Km-\)P-\)'\^ ^ first column, and 

x{Np-\)^^k,m^\^^u{Np-\)^'"^^k,m^M~\^^u{Np'\)\ its fifst TOW, whetc {c/^^,} ^Te thc 
binary phase shift keying (BPSK) data bits of the wireless terminal. 

In addition, the temporal signature matrix Sf^\L'K{N ^-Q-\)p is defined as 

10 = G^Sjt . The row of Sj^ denotes the temporal signature of the /^^ path of the h!^ 
wireless terminal seen at the receiver. Similarly, the colunm of Ak denotes the spatial 
signature of the path of the 1^ wireless terminal. If we consider the impulse response 
of the wireless channel, and assume there is only one symbol transmitted for each 
wireless terminal, then Dk equals dk^\ times an identity matrix. After stacking the 

1 5 columns of X, equation (2) can be written as 

vec(X) = Y.d,M^K)Pk (3) 
= Z^*.. Z >^*.'^(^*,. ® ^*./) + ^ec(N) (4) 

where the "vec" operator stacks colunms of a matrix putting them one below another, ® 
denotes the Kronecker product, o denotes the Khatri-Rao product (column-wise 

20 Kronecker product), and Sj^ / , is the row of 5^ . Thus, S^j ® a(^j^ /) is the space-time 

signature of the path of the //^ wireless terminal. If there is more than one symbol 
transmitted for each wireless terminal, inter-symbol interference results in space-time 
signatures overlapped in time v^thin the received signal X. 

Returning now to Fig. 1, the joint channel estimator 1 14 receives the baseband 
25 signal and, using the above described model, determines characteristics of the channels of 
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the wireless terminals 110. Channel information including the path DO As (direction of 
arrival) { O^j}^ propagation delays { ^ }, and the complex path fading amplitudes 
{ fii^ j } are required in the space time matched filter 1 16, the signal detector 120, and the 

code optimizer 118. The joint charmel estimator 1 14 determines these channel 
characteristics as follows. 

The joint channel estimator 1 14 estimates the path DO As (^^^ ,) and the 

propagation delays (r^,) by using known subspace-based algorithms. Appropriate 

subspace-based algorithms which may be used by the joint channel estimator 114 are 
described in, M. C. Vanderveen, C. B. Papadias, and A. Paulraj, Joint Angle and Delay 
Estimation (JADE) for Multi-Path Signals Arriving at an Antenna Array, IEEE 
Communications Letters, vol. 1, no. 1, pp. 12-14, Jan. 1997. Since both the path DO As 
and the propagation delays usually remain unchanged over a long period of time, the 
above referenced subspace-based algorithms may be used in order to obtain accurate 
channel estimations. For purposes of this description, it is assumed that we have perfect 
information on the path DOAs and the propagation delays as a result of one of the above 
referenced algorithms. With the path DOAs and propagation delays available, { Aj^ } can 
be reconstructed by using the array manifold a(.) and {Gk} can be reconstructed by using 
the temporal manifold Note that {Sk} denotes the wireless terminal codes, which are 
known to the base station 102. The only remaining unknowns left in equation (2) are 
{Bjt} and {Djt}, which are the complex path fading amplitudes and the transmitted data 
symbol matrix respectively. The data symbol matrix can be obtained by using known 
training sequences or in a decision-directed fashion, i.e., the demodulated data symbols 
can be easily used as training data. To estimate {Ek} , we can re-write equation (2) as 
X=[Ai, A2, ...,A/^]diag(Bi,B2, Bic) diag(Gi, G2, ...,Gk) 
diag(Si, S2, Sj,) [D r 1 , D D ] ^-h N (5) 
= ABGSD + N, (6) 
where A = [Ai, A2, Aa:]; B = diag(Bi, B2, Bj^); G = diag(Gi, G2, Ga:); 
S = diag(Si, S2, . . . , Sk); and D = [D[1,D[,...,D^]^. The least square criterion is used 
to estimate the fading amplitudes hat as 
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P = ((GSD)^ 0 A) \ec(X) (7) 

where P = [^[>^[,".,^j^]^ with being the least square estimate of ; and (.) ^ 
denotes the pseudo-inverse. 

The channel estimates calculated by the joint channel estimator 1 14 are provided 
5 to space-time matched filter 1 16, code optimizer 118, and signal detector 120. 

The space-time matched filter 1 16 is used to obtain sufficient statistics for 
demodulating the digitally modulated signals. The output of the space-time matched 
filter 1 16 may be written as 

X,,„ = Re{/r[B;A;X'"S;}}. (8) 

10 Where "tr" denotes the trace of a matrix, and X" : P x (N + Q - V)p is copied from the 

({m-iyp+l)'" to {mNp + Qp) columns of X. The matrix multiplication of A\ and 

X'" denotes the path- wise beam-forming. Let X^ ^ = A*^X'" . The trace of the matrix 

multiplication tr { B^Xj^^Sj^ } denotes the matched filter for the temporal signature and 

the maximum-ratio path combining of the wireless terminal. To expand X^^^ 
1 5 equation (8) can be rewritten as 

/=/ 

where N'" and D " are the sub-matrices of N and Dj respectively, so X'" is constructed 
fi:omX. Note that is a Toeplitz square matrix with [^,>^Oi^(^^j),£/, „,_, Oi^^(^_,^^_j j]^ 

as its first column, and [<2f,^;„^0i^(jv^,),£/, O,^^^^.,^^, as its first row. Therefore, the 
20 first term on the right hand side of equation (9) can be written as 
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= t Rew,.„>9;((A;A,) o (s,s;)0^, +c/,,„.X((a;a,) 0(s,(-hi)s;(-i))'oa 
+ d,^^X(iK^>) o (s,(-i)s;(+i))^);9,}, (10) 

where O denotes the Hadamard product or the entry- wise product, S4(+l) and (-1) 
are both sub-matrices of . The rightmost and the left most {Q-\) p columns are 

respectively denoted by Sji(+1) and Sj(-l). For additional information on the 
Hadamard product/entry-wise product, see Roger A. Hum and Charles R. Johnson, 
Topics in Matrix Analysis, Cambridge University, Press 1991. 

We can stack the matrices and write the received signal at the output of the space- 
time matched filter 1 16 as 

X= i'd + N, 

(11) 

where X = [X, , Xj,, , . . . , , X, 2 , . . M Xj^. AY ]^ d = , ^2 p . . . ,£/jc p £/^^^ 

N = [N,,,,N2.p...,N^,N,2,...,N^,;^,]'' with N,^=Re{/r{B;A;N'"S:}}. If we define 
H(0): 

/^xA; H(+l):/:x/:, andH(-l):/Cx Kbs 

H(0) = B*((A*A) O (SS*)^)B (12) 

H(+1) = B'((A*A)0 (S(+1)S*(-1))^)B, (13) 
H;-1)=B*((A*A)o (S(-1)S*(+1))'')B (14) 

with 

A = [A, Aj,..., A^], S = [S^SJ,...,S^]^S(+l) = [Sf (+1),S[(+1),...,S^ (+l)]^ 
S(-l)=[Sf(-l),S[(-l),...,S^ )(-!)]" 
and define the fading amplitude matrix B:KLy.K as B = diag {{P^pj...,^,.) , then 
can be expressed as 
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Re{ 



e(0) H(-i) 

H(+l) e(0) H(-l) 
H(+l) H(0) 



H(-l) 



0 



} . 



(15) 



0 



e(+i) e(0) 



The output of the space-time matched filter 1 16 is provided to the signal detector 
120. The signal detector 120 also receives the channel estimations from the joint channel 
estimator 1 14 and the codes which are output from the code optimizer 118. We assume 
that only in-cell wireless terminals need to be demodulated at the base station 102. 
Therefore, only the space time matched filter 1 16 outputs related to these ATo wireless 
terminals are processed in the signal detector 120. We now describe three possible 
embodiments of the signal detector 120. 

First, signal detector 120 may be implemented as a single wireless terminal 
detector. In this embodiment slicers are used at each output of the space-time matched 
filter 116. Decisions are instantaneous with no delay. The data estimate is given by 



where X^ ^^ is the output of the space time matched filter 116 in equation (8). 

In a second embodiment, signal detector 120 may be implemented as a decorrelator. In 
this embodiment, a zero-forcing solution is applied to solve equation (11). The data 
estimate can then be written as 



In a third embodiment, signal detector 120 may be implemented as a minimum mean 
square error (MMSE) detector. This embodiment generally performs better than the zero- 
forcing solution of the decorrelator in the second embodiment, when taking into account 
the effect of the correlated noise. The coloring matrix of the noise, N in equation (11), 
at the output of the space time matched filter 1 16 is given by: 




sgn(X,^J, 



(16) 



(17) 




% 



(18) 
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The data estimate for the MMSE detector can then be expressed as 



(19) 



where al denotes the noise variance chosen to match the desired — with £^ being the 

bit energy and A^^ being the average noise variance, and cr^ = 1 for the BPSK 
modulation scheme. It is noted that in both the decorrelator and the MMSE detector 
embodiments, the maximum possible delay could be up to A/ symbols because of the 
block processing of the asynchronous signals. The reverse link data outputs of the signal 
detector 120 may be provided to, for example, a telephone network node such as a 
switch. 

In accordance with the invention, the base station 102 includes a code optimizer 
1 1 8 for assigning codes to wireless terminals. The steps performed by the code optimizer 
1 1 8 will be described in conjunction with the flowchart of Fig. 2. In one embodiment, 
the code optimizer 1 18 is implemented using a computer processor executing stored 
computer program code. Thus, the code optimizer operates under control of the 
processor executing stored computer program code. The computer program code defines 
the steps described herein in conjunction with the flowchart of Fig. 2. Given the 
description contained herein, one skilled in the art could readily implement the code 
optimizer 1 18 in other ways. For example, the code optimizer 118 could also be 
implemented using hardware, or some combination of hardware and software. 

The steps shown in the flowchart of Fig. 2 are performed by the code optimizer 
1 18 in order to optimize the codes assigned to the wireless terminals communicating with 
the base station 102. In step 202, as processing set is initialized with one of the wireless 
terminals. In step 204 a target wireless terminal is chosen from the processing set. If 
step 204 is executed immediately following step 202, then there will only be one wireless 
terminal in the processing set and that wireless terminal will be chosen as the target 
wireless terminal. In step 206 a random code search is performed for the wireless 
terminal chosen as the target in step 204. The random code search is performed in order 
to try to improve the current code assigned to the target wireless terminal by attempting 
to reduce the interference experienced by the target wireless terminal. 
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One way to suppress the mutual interference among the wireless terminals is to 
maximize the minimum signal to interference-plus-noise ratio (SINK) among in-cell 
wireless terminals, where 

5 The SINR of the k!^ wireless terminal can be written as the power ratio of the signal 
and the interference-plus-noise: 

^*=y^fV- (21) 

The signal power <^^:Kx \ can be written as 

= diag(M{0)) O con/(diag(H(0)), (22) 

10 where H(0) is given in equation (12), "conj" is the complex conjugate operation, and the 
output of diag(.) is a colunm vector consisting of the diagonal elements of the argument 
matrix. The interference power <^j:Kxl can be written as 

4,=diag{im')-<^s^ (23) 
where H = [H(-1),H(0),H(+1)], and H(0), M (+1), and H[(-l) are given in equations 

15 (12), (13) and (14) respectively. The noise power (^j^:Kxl can be written as 

4,=^diagmO)y (24) 

Note that in using the minimum SINR as the criteria for code optimization, two factors 
are taken into consideration: 1) the received power of each wireless terminal at the base 
station, and 2) the orthogonality among codes. 

20 Although some out-of-cell wireless terminals may generate interference with 

respect to in-cell wireless terminals we assume that the only visible wireless terminals are 
in-cell vmeless terminals such that K = K^, For systems with a high frequency re-use 
factor this is a reasonable assumption. 

We now retum to step 206 of Fig. 2 in which the code optimizer 118 performs a 

25 random code search. The less dense regions of the signal space must be explored in order 
to maximize the distances between the wireless terminal codes. More particularly, we 
need to adopt some of the un-used codes by random search in order to span the entire 
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1 



signal space. It is also noted that a good code for one wireless terminal may not 
necessarily be a good code for another wireless terminal because the wireless channel of 
each wireless terminal is different. Therefore, a random-searching procedure of codes is 
required for each wireless terminal in order to cover those less dense areas in the signal 
5 space. When the number of wireless terminals is small compared to the dimension of the 
signal space, it is easier to find a less dense area by random searching. Otherwise, longer 
searching may be required. We now describe two alternate techniques for performing the 
random code search step 206. The first technique is to continue searching all the possible 
codes until a better code is obtained. If the number of wireless terminals is large 
10 compared to the dimension of the signal space, it may take a long time to find a better 
code. The second technique is to only search a subset of randomly-picked codes, and 
choose the best code in the subset if it is better than the original code prior to the random 
2 search. If the number of wireless terminals is small compared to the dimension of the 

■=y signal space, the second approach might require a large amount of urmecessary searching. 

iil 

\d 1 5 However, if the nimiber of wireless terminals is large compared to the dimension of the 
;E signal space, this second technique might stop at a sub-optimal value. One skilled in the 

□ art would recognize that other techniques could be used to implement the random code 

■=^1 search. 

Although it is unlikely that the best code will be found by the random search step 
\M 20 206, the best available code is likely to be in the vicinity of the less dense area in the 
, J signal space which surrounds the code found in step 206. Thus, the code for the target 

wireless terminal may be further optimized by gradient-searching this less dense area in 

the signal space. 

In step 208 it is determined whether the random code search of step 206 resulted 
25 in an improvement of the SINR for the target wireless terminal. If there was an 

improvement, then in steps 210 and 220 two types of gradient searches are performed in 
order to find the best code for this wireless terminal. 

In step 210 a gradient search of codes is performed as follows. A small change is 
made to the code of the target wireless terminal to determine if there is improvement in 
30 terms of the minimum SINR. Assume the length of the code is N chips long. A'^ different 
codes are generated by individually changing the sign of each chip in the code. These 
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different codes are then compared with the original code. If none of the new codes 
performs better than the original code, then the gradient search of codes is considered to 
be converged. Otherwise, the original code is replaced with the best code among the 
new codes, and the same process is repeated. Note that the change of the sign of only one 
5 chip in the code results in only a minimum change of the code on the transmission side. 
It is expected that such a minimum change in direction will also result in only a minor 
change of the code as seen at the receiver. 

In step 220 a gradient search of transmission delays is performed as follows. The 
delay of transmission of the wireless terminal code is adjusted by"a small step AT" in 

10 order to determine the optimum delay of a specific code. Assuming band-limited signals, 
the step size AT" should be in inverse proportion to the signal bandwidth. 

As represented by step 222, the gradient search of codes (step 210) and the 
gradient search of delays (step 220) continues until no further improvement in the SINR 
for the target v^reless terminal is achieved. It is determined in step 224 whether the code 

1 5 searching steps (steps 204-222) have converged. Convergence of the code search is 
defined as the condition where there is no space for further improvement regardless of 
which target Mdreless terminal of the processing subset is picked. Convergence may be 
determined by repeating steps 204 through 208 until there is no improvement. Eventual 
convergence is guaranteed because there are only a finite number of possible code 

20 assignments, although this finite number 2^G^ could be very large, where G is defined 
as the number of possible transmission delays. There are a finite number of possible 
transmission delays because of finite resolution in time for the assumed band-limited 
signals in the system. If an exhaustive search is possible, the optimum assignment is 
guaranteed. Otherwise, since we are constantly looking for a better minimum SINR, the 

25 proposed algorithm can at least converge a local optimal solution. 

Returning to Fig. 2, if the test in step 224 is no, then control returns to step 204 
and a new wireless terminal firom the processing set is chosen as the target wireless 
terminal and the code search of steps 206-222 are repeated. Note that the new target may 
be a wireless terminal already chosen during this processing loop. If the test in step 224 

30 is yes, then in step 226 it is determined whether the number of wireless terminals in the 
processing set equals the total number of wireless terminals conmiunicating with base 
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station 102. If not, then in step 228 another wireless terminal is added to the processing 
set and steps 204 through 224 are repeated for the new processing set. If it is determined 
in step 226 that the number of wireless terminals in the processing set equals the total 
number of wireless terminals communicating with base station 102, then in step 230 the 
codes which are currently assigned to the wireless terminals are output from the code 
optimizer 118 and processing ends. 

Retuming now to Fig. 1, the optimized codes from code optimizer 1 18 are 
provided to the signal detector 120 and to the forward link modulator 122. The forward 
link modulator 122 also receives forward link data inputs which may be, for example, 
telephone signals from a telephone network. The forward link modulator 122 modulates 
the forward link data and the codes received from the code optimizer 1 1 8 and provides 
the modulated signal to the frequency up-converter 124. The frequency up-converter 124 
transmits the signal to the wireless terminals 1 10 via antenna 126. In this manner, the 
newly assigned codes are transmitted to the wireless terminals 1 10. The wireless 
terminals then use the newly assigned codes for subsequent reverse link communication. 
As described above the signal detector 120 also receives the new codes from the code 
optimizer 1 18 so that the signal detector 120 can decode the subsequent reverse link 
conmiunication. 

Thus, as described above in connection with one advantageous embodiment, the 
use of a processing set results in an iterative algorithm in which wireless terminals are 
added into the processing subset one by one and only after the code optimization 
procedure has converged for the prior processing set. This technique, which we call 
sequential packing, is advantageous because the code assignments are more likely to 
converge to the optimum solution when the number of wireless terminals is small. 
Furthermore, the optimal code assignment of a A^-wireless terminal system is likely to be 
similar to that of a (A^-l) v^reless terminal system. 

Further, during the processing of each processing set, we use the technique of 
target wireless terminal isolation. Instead of updating all the codes for all the wireless 
terminals in the processing set, we isolate a target wireless terminal each time and 
optimize its code before moving on to optimize other wireless terminal codes. This also 
speeds up convergence because before the convergence of the code optimization, most of 
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the codes may have already achieved good orthogonal distribution. Disturbing all the 
codes would likely slow down the convergence speed. 

The following portion of the description addresses performance issues of the 
above described code assignment techniques. The performance of the algorithm 
5 described in conjunction with Fig. 2 depends on the number of wireless terminals, the 
dimension of the signal space, the overall channel of each wireless terminal, and the 
adaptation procedure used for code optimization. 

In general, there are 2^G^ possible code assignments, where G is defined as the 
number of possible transmission delays, K is the number of wdreless terminals, and is 
10 the length of a code. There are a finite number of possible transmission delays because 
of the finite resolution in time for the assumed band-limited signals in the system. The 
minimum SINR can be calculated for each code assignment. System performance is 
O improved in accordance with the described technique of replacing the original random 

ilj signature assignment with one of the best signature assignments if possible, and thus 

^fl 15 improve the minimum SE^R. Fig. 3 shows an example of the minimum SINR 
£ distribution in a simple system in which we assume a synchronous system with only one 

ill 

path of unit gain for each wireless terminal, i.e., G = 1 , a code of 4 chips long, up to 4 

L E 

Q wireless terminals, a single antenna, and — = oo. Each curve in Fig. 3 is labeled by the 

i7j number of wireless terminals in the system. J is defined as the total number of possible 

'5 20 code assignments, and therefore, J = 256, 4096, and 65536, for two, three, and four 

wireless terminals respectively. In Fig. 3 the minimum SINR vs. x is plotted, where the 

minimum SINR is sorted in ascending order, x = — , and J is the code assignment number 

after sorting. Therefore, the horizontal axis of Fig, 3 represents the percentage of the 
bottom performance among all the possible code assignments. Note that the minimum 
25 SINR is very low most of the time without code assignment optimization regardless of 
the number of wireless terminals in the system. The improvement on the minimum SINR 
caused by reassigning the wireless terminal codes is very significant when the number of 
wireless terminals is close to the dimension of the signal space. When the number of 
wireless terminals is small, the improvement is still great but not as significant, because it 
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is relatively easy to separate the code vectors in the signal space. In Fig. 4, the top 2 
percent performance of Fig. 3 is expanded. In this example, we can obtain orthogonal 
codes if the algorithm of Fig. 2 converges to the minimum SINR among the top 20 
percent, top 1 .5 percent and top 0.05 percent, respectively, when there are 2, 3 and 4 
5 wireless terminals in the system. 

The problem of how many searches are required to reach a specific percentage of 
the top minimum SINR performances is now addressed. The trade-off between the 
complexity of the search algorithm and the improvement on the SINR performance is 
also investigated. Since the number of possible code assignments 7 is usually a very 
10 large number, the efficiency of the code searching turns out to be an important issue. 

First, we analyze the probability of each code assignment getting picked using two basic 
strategies: 1) blindly searching all the possible codes assignments until a signature 

0 assignment with better minimum SINR is found, and 2) blindly searching J\ code 
y assignments and picking the best one. 

J";-: 15 We estimate the approximate number of blind searches in the algorithm and then 

P compare this number to the number of blind searches required for the minimum SINR to 

ill 

Q reach a certain performance level. We first analyze the code assignment searching 

strategy in which all the possible code assignments are blindly searched until an 
assignment which achieves a better minimum SINR than the current assignment is found. 

1 jj 20 /7Q(y) is defined as the probability of picking the signature assignment. For a system 
% using pseudo random codes, the probability distribution is given as 

Poij) -—J - 1,2,..., J. Recall that the code assignments are sorted in ascending order. 
J 

If the probability that the y>l assignment gets picked next is zero, then the probability 
that any of the last J -71 assignments gets picked next is — ^- — . Therefore, the 

25 probability distribution after the first cycle of searches becomes 

p, =Mpo, (25) 

where Po =bo(A/?o(-^-l) AO)]^Pl =[A(AA('^-l)>-.A(l)r' and the 

transition matrix M : / x J is defined as 
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10 



15 



20 



1 1 i i 

2 3 

0 0 i i 

2 3 

0 0 0 - 
3 



0 0 0 0 
0 0 0 0 



1 



J-\ 
1 

J-1 
1 

J-1 
1 

J-1 

0 



(26) 



From observing the structure of M, it is seen that this code assignment searching 
algorithm pushes the distribution of the minimum SINR toward the better side of 
performance. To further improve the minimum SINR, we can repeat the same searching 
strategy several times. If / search cycles are conducted, we can obtain a probability 
function p, given by 

p,=M'Po (27) 
Note that once a good minimum SINR has been obtained, the number of searches 
required in each search cycle increases as the searching continues, since it is getting more 
difficult to find an even slightly better minimum SINR. The average nimiber of blind 
searches required at the / cycle can be written as 



"/=ZaO)4-^+(^-i)a(</) 



(28) 



The number of possible signature assignments J is usually a large number. Therefore, 
equations (25), (27) and (28) are difficult to evaluate. Thus we need to rewrite these 
equations in a continuous form. Let the assignment be the bottom x(j) percent among J 

assignments, namely x = — , and let p.{x) be the probability density function of x after 

the search cycle. The probability density function after the first cycle of searching 
becomes 



1 



-dx = In- 



1 



1-x {\-x) 

where "In" denotes the natural log. Similarly, the probability density function after the 
cycle of searching becomes 



(29) 
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m 



PM = llp,.,(x)-^dx. (30) 
\-x 

By solving equation (30), we obtain 

Mx) = ^(ln-^y, (31) 

where /! denotes factorial i. The cumulative distribution function (CDF) of the minimum 
5 SINR after the /'^ cycle of searching can then be written as 

A .CDF(x) = X - 2;-bil(/«-L)* . (32) 

In addition, the average number of blind searches in the i'^ search cycle can be written as 

=loA-iW7^^- (33) 

^ Fig. 5 shows the improvement of the probability density function of the minimum 

g 10 SINR using the random search technique of searching all possible codes until a code 

which is better than the current code is found. Each curve in the figure is labeled by the 
number of the search cycles. By comparing Fig. 3 with Fig. 5, the number of cycles of 
i ji blind searches which are required to achieve a certain minimum SINR performance can 

be determined. Fig. 6 shows the number of blind search cycles required to reach certain 
3 15 performance levels for the two random search techniques. Section 602 of Fig. 6 shows 
^ the number of search cycles required to have 50, 90, and 99 percent chance for the search 

^ results to be among the top 10, 1 and 0.1 percent of performance of all possible code 

•0 assignments for the first random search technique which searches until a better code is 

found. For example, to have a 90 percent chance to reach the top 1 percent of 
20 performance using the first random search technique, we need to search for a better 
assignment at least 7 times. 

As described above, the second search strategy is to randomly search through a 
fixed number of code assignments to find the code with the best minimum SINR. If we 
search through / possible code assignments, the probability of picking the assignment 
25 PiU) is given by 



■for j>i 



Ci . (34) 

0 for j < i 
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However, if we use i - 1 pairwise comparisons to search through i possible code 
assignments and we allow the same code assignments to repeat within i assignments, we 
can modify p^ij) as follows: In each pairwise comparison, if the assignment is first 

used, there is a ~ — ^- chance that the assignment will be compared to an assignment 



J-\ 



with poorer minimum SINR, while there is a 



1 



J-1 



chance that the assignment will 



be compared to an assignment with better minimum SINR. Therefore, the probability of 
the first y'l - 1 assignments getting picked next is zero; the probability of the same 

assignment getting picked next is -^-^ ; and the probability of any of the better 

*J 1 



assignments getting picked next is 



1 



J-1 

'-I, 



Thus, we can write {/?, } as 



10 p,=M'-'po.' = 2,3,... (35) 

where Po =[a(Aa(-^-1X-,/?o(0]^P, =[a(Aa(-^-1X-,;7,(1)]^ and the 



transition matrix M : J x J is defined as 

1 1 



M= 



1 



1 



J-\ J-1 
J-2 1 



J-1 J-1 
1 1 



J-1 J-1 
J-3 



0 0 

0 0 
0 0 



J-1 

0 
0 



J-1 J-1 
1 1 
J-1 J-1 



1 



1 



J-1 J-1 

0 0 



(36) 



1 5 Again, the structure of the matrix M shows that this code assignment searching algorithm 
pushes the distribution of the minimum SINR toward the better side of the performance. 
Because J is usually a large number, equation (35) is not easy to evaluate. Therefore, pi 
is expressed in a continuous form as 

Pi{x) = \lpQ{x)dx + p^ix)x = lx (37) 
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assuming Pq{x) = l,jc €[0,1], and also 



Pi (^) - io A-i (x)dx + Pi,, / = 3,4. . . 



(38) 



By solving equation (38), we obtain 



(39) 



The CDF of the minimum SINR after searching through i signature assignments can then 
be written as 



Fig. 7 shows the improvement of the probability density function of the minimum 
SINR using the second search strategy of searching through a fixed number of code 
assignments. Each curve in the figure is labeled by the number of assignments searched. 
Comparing Fig. 3 with Fig. 7, we can determine how many code assignments must be 
searched in order to achieve a certain minimum SINK performance. Section 604 of Fig. 
6 shows the number of searches required to have 50, 90, and 99 percent chance for the 
search results to rank among top 10, 1 and 0.1 percent of performance of all possible code 
assignments using the second random search strategy. 

There is a trade-off between the minimum SINR performance and the 
convergence speed. In the described algorithm, the convergence criterion, the step-size 
of the gradient search, and the technique of picking a target wireless terminal can be 
controlled in order to adjust the algorithm complexity and convergence speed. If the 
complexity of the assignment searching procedure is increased, then the minimum SINR 
performance is expected to improve, while it might take more time for the algorithm to 
converge. This trade-off is investigated in order to determine how much increase in 
algorithm complexity is required to achieve a certain performance improvement. The 
number of equivalent blind searches h in the proposed algorithm is estimated. Based on 
fi , and Figs. 5-7, the improvement on the minimum SINR among all the possible code 
assignments is evaluated. 

As described above, blind searches are generally inefficient regardless of the 
searching strategy used. In view of the above performance analysis, improvements for the 
searching techniques may be summarized as follows. 



p,,CDF(x) = x'. 



(40) 
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The random searches and the gradient searches serve the respective purposes of 
region dividing and local regional searching in the signal space. The random search first 
looks for the regions where a higher minimum SINR is more likely to occur. This 
random search may be seen as dividing the signal space into small regions and then 
5 finding one of the best regions among other sub-optimal regions. The gradient search 
then searches for the best minimum SINR around the best local region determined by the 
random search. The number of random searches for the code of wireless terminal k is 
defined as h^ j^ . If the best code assignment in the best region has higher minimum SINR 

than any code assignments within those sub-optimal regions, then the equivalent number 

1 0 of blind searches for wireless terminal k can be written as 

«.=^.A>(A' + 2). (41) 
^ Equivalent gradient searches are assumed to be conducted for each of the sub-optimal 

.p regions although these gradient searches are actually omitted. Note the factor A'^ + 2 in 

% equation (41) represents the number of searches in each step of the gradient search; A'' 

1 5 searches for the code, and 2 searches for the transmission delay, 
ji With respect to target wireless terminal isolation and sequential packing, if the 

^ number of wireless terminals in the system is much less than the dimension of the signal 

□ space, it can be assumed that the code search for each wireless terminal is independent. 

^ In other words, we assume that the best code found while searching the best codes for 

=2 20 wireless terminal k regardless of the codes of other wireless terminal, is also the best code 

u 

Q for wireless terminal k when searching the best code assignments for all the wireless 

terminals. The main purpose of the sequential packing technique described above is to 
validate this assumption. Therefore, the overall equivalent number of blind searches h in 
the algorithm is roughly given by 

25 n = f[n,=f\h,,h^,{N + 2) (42) 

The estimation of fi in equation (42) is over-optimistic; however, it offers some insights 
into the relationship between the system performance and the algorithm complexity. 

The foregoing Detailed Description is to be understood as being in every respect 
illustrative and exemplary, but not restrictive, and the scope of the invention disclosed 
30 herein is not to be determined from the Detailed Description, but rather from the claims 
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as interpreted according to the full breadth permitted by the patent laws. It is to be 
understood that the embodiments shown and described herein are only illustrative of the 
principles of the present invention and that various modifications may be implemented by 
those skilled in the art without departing from the scope and spirit of the invention. 
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